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Abstract 

We report the first lattice QCD results of the scattering amplitudes of the Kir system for I = 1/2 
channel together with 7 = 3/2 case. We investigate all quark diagrams contributing to these iso- 
spin states, and find that the scattering amplitudes are expressed as combinations of only three 
diagrams after setting the masses of u-quark and ci-quark to be the same. 

The lattice simulations are performed in the quenched approximation at (5 = 2.23 on a 12 3 x 24 
lattice with an improved Iwasaki gauge action. We employ a new dilution-type noise method to 
get accuracy of data with reasonable CPU time. A simple method is proposed and applied to 
eliminate lattice artefact due to the finite extent of lattice along the time direction. 

A clear difference in the quark mass dependence between 1 = 3/2 and I = 1/2 channels is 
observed. Although the chiral extrapolation is subtle, we assume E] <n oc vr? u d , and obtain the 
5-wave scattering lengths as a (I = 3/2)m n = -0.084^;^ and a (I = l/2)m 7r = -0.625 ± 0.012. 
We show all necessary formulas which make the calculation possible. 

We argue that AN is the most appropriate target of the Liischer's formula for baryonic system 
because it has no it exchange diagrams and has a scattering length suitable for a lattice QCD 
simulation. 

PACS numbers: 12.38.Gc, 13.75.-n, 21.30.Fe, 24.85. +p 



*Present address: Admissions Center, Hiroshima University, Higashi-Hiroshima 739-8511, Japan 
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I. INTRODUCTION 



One of the main objectives in our study of hadron nuclear physics is to describe hadron 
interactions on the basis of QCD (quantum chromodynamics). The study of hadron inter- 
actions will help us understand such interactions in terms of multi-quark reactions mediated 
by gluons. This is still simply a dream, and we need using phenomenological models whose 
parameters are fitted to the experimental data. However, even this strategy is still diffi- 
cult, particularly when the strangeness is included, because in such a case, there is limited 
experimental information available. 

The development of high-energy hadron accelerators such as those at JLab, LEPS, and J- 
PARC and hadron studies carried out at these facilities have contributed to the accumulation 
of experimental data pertaining to high-energy quark interactions. In these accelerators, not 
only u and d quarks but also s quarks are excited, and high-statistics studies on hadronic 
reactions have been carried out; numerous hypernuclei have also been produced in these 
reactions 

The study of hadronic reactions in a unified manner through first-principle calculation is 
very important, then, quark-gluon reactions can be studied on the basis of QCD by analyzing 
the large amount of precise data obtained at the aforementioned facilities. The results of 



recent lattice calculations of the NN force are very encouraging in this direct ion [lj. We will 
step into a new era of hadron physics. 

Liischer derived the basic formula for the calculation of scattering lengths on the basis 
of lattice QCD simulations, where the S"-wave scattering length between two hadrons is 
related to the energy shift of the two-hadron state that is confined in a finite periodic spatial 
box of size L 3 at zero relative momentum Meson-meson, meson-baryon, and 

baryon-baryon scattering lengths have been studied recently by using Luscher's formula f], 



3, lSl, la, llO, llll, • In these calculations, contributions from different types of quark diagrams 
to the hadron four-point functions have been analyzed. Luscher's formula is expected to play 
a major role in the study of hadronic reactions on the basis of lattice QCD. In principle, 
there are no limitations for including strangeness and other flavor degrees of freedom, and 
no free parameters. 



We have calculated scattering lengths in the Ap system |9( using Luscher's formula. This 
method involves a long calculation time. To reduce the calculation time, we use the modified 




n 
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noise method in the present study for the evaluation of quark propagators. 

In this research, we study Kir scattering, which is a simple but important fundamental 
reaction. This reaction is interesting because of reasons as follows : 

1. It is the simplest reaction that includes s quark. 

2. The 1=1/2 channel of this reaction is directly related to the scalar meson k [3]. 
Further, this reaction is easier to study than the irn scattering reaction, in which case 
I = is directly related to the scalar meson a. 

3. The force between K and tc may produce a K-kN bound state, which can be used to 



explain the penta-quark state 14], ll5j 



4. Direct lattice QCD measurement of I = 1/2 and 3/2 will provide a test of the validity 
of chiral perturbation with strangeness. 

In the present study, we have evaluated the scattering length of the Kir system by lattice 
QCD owing to the above mentioned features of Kir system. 

In Section II, we explain the formulation based on Liischer's formula for the Kir system. 
In the K7r system, there are 22 quark diagrams. The number of independent diagrams can 
be reduced to six if we assume the masses of u and d quarks to be identical. After simple 
calculation related to the iso-spin states, only three diagrams contribute to the 7 = 1/2 state 
and two to the 1 = 3/2 state. In Section III, we show the results obtained in our simulations 
for the iso-spin channels 1 = 3/2 and 1/2. In Section IV, we discuss the differences in the 
contributions from each diagram. The final section includes concluding remarks. 

II. METHOD 

A. Scattering length determined using Liisher's formula 

Liischer's formula which relate the energy shift AE to the scattering length|2| is given as 
AE = E Kn - (m K + m*) 

_ 2^^ [1+ei?+C2(g)2]+0(L _ 6) (1) 

with ci = —2.837297 and C2 = 6.375183, where Ek-k is the total energy of Kit system, uik 
and are masses of K and ir, and L represents the spatial size of lattice, respectively. 
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Rummukainen and Gottlieb extended the above formula to moving frames [16( and suc- 
ceeded in calculating phase shifts in addition to the scattering length. 

Using operators Ok(x\) and O n (x 2 ) for K and tc at points X\ and x 2 , respectively, we 
represent hadron four-point functions as follows : 

C Kn (x' 1 ,x' 2 ,x 1 ,x 2 ) = (OKix'jO^x'jOlcixjOKxz)). (2) 

Here, (• • •) represents the expectation value of the path integral, which we evaluate using 
quenched lattice QCD simulations. 

After obtaining the sum over spatial coordinates x%, x 2 , x[ and x 2 , we obtain the four- 
point function in the zero-momentum state, whose behavior is given below : 

Y Y Y Y C^K, x' 2 , x h x 2 ) 

33i 3*1 3*2 

= Z Kir cosh(E Kw (t - N t /2)) + Z' Kv cosh(^(t - N t /2)) + • • • . (3) 

Here, x[ = (x[,t' 1 ), x' 2 = (x 2 ,t' 2 ), x\ = (xx,ti), and x 2 = (x 2 ,t 2 ) with t[ = t' 2 and t\ = t 2 . t 
stands for the time difference, t = t 2 — t\. E KlT and E' Kir are the ground and excited levels, 
respectively. Hadron two-point functions are also given by 

Ck{xx) = (Y°k( x 'i)Y°k^) 

x[ xi 

= Z K cosh(m x (t - N t /2)) + Z' K cosh(m' K (t - N t /2)) + ■■■ 
C„{x x ) = (Y0^2)Y ^) 

= Z^ cosh(m^ - N t /2)) + Z' w cosh{m' n {t - N t /2)) + ■■■ . (4) 

The energy shift AE can be deduced directly as the difference between Ex n and Ex + E n , 
which are obtained from our simulations of the four-point function Ck-k and the two-point 
functions, Ck and C n . In the present study, the final values of AE are obtained from the 
results of the fitting procedure for Ck-k, Ck and C n . The details will be explained in Section 
III. 
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B. Quarks diagrams in Kit scattering 



In the Kit system, there are two iso-spin states, 1 = 3/2 and / = 1/2, 

\Kit{I = 1/2) > = \ \\K° > |tt + > — ]=\K+ > |tt° >, (5) 
V 3 \/3 

\Kit(I = 3/2) > = \K + > |tt+ > . (6) 
Then, the scattering amplitude M for Kit scattering in the 1 = 1/2 state is given by 

M(I = 1/2) = | < K°7T + \S\K°7T + > < K°7T + \S\K + 7T° > 

_ yl < K + 7T°\S\K°7T + >+-< K + 7T°\S\K + 7T° > (7) 

3 3 
We introduce the following operators explicitly for Ok(x) and O n (x) : 



K o(x) = s(x)j 5 d(x), 

K +{x) = s(x)"f 5 u(x), 

On+(x) = — J(x)7 5 m(x), 

0^0(2;) = -j=^u{x)^u{x) - d(x)~/ 5 d(x)y (8) 

We insert Eq. (8) into Eq. (2) to obtain the quark diagrams for Kir scattering. The details 
of our calculations are provided in Appendix [K\ in Eqs. Al - A22. Finally, we have 22 
different diagrams for / = 1/2, as shown in Figs. [Hand [2] , which correspond to Eqs. Al - 
A22. 

However, by assuming that u and d quarks have the same mass, we can categorize the 22 
diagrams into six independent groups. Diagrams 1, 12, and 22 in Figs. [Hand [2] are compiled 
into Group 1. Similarly, diagrams 11, 17, 19 and 21 are compiled into Group 2, No. 2, 4, 8, 
and 13 in Group 3, No. 6, 10, and 15 in Group 4, diagrams 3, 5, 14, and 18 to Group 5, and 
diagrams 7, 9, 16, and 20 Group 6, respectively. According to Eq. (8) and Eqs. Al - A22, 
the weights of each of these groups are given as follows: 



Group 1: 

I + + = 1:4 

3 3 V 2 A ' 3 V 2 A ; 

Group 2: 
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FIG. 1: Quark propagators corresponding to Eq. in Appendix A. 



Group 3: 
Group 4: 
Group 5: 
Group 6: 



(9) 



|)(-i) + (-f )(--!) + <-f )(-^)(-i) + = -§ = *, 

-^)(-^)( +1 ) + (-f)(-L)(- 1) + (i)(-i) = i:X, 



Here, A, if, and X are denoting the type of the quark diagrams shown in Fig. [31 Finally, only 
the three diagrams shown in Fig. [3] remain to contribute to the Kn scattering amplitudes. 
Ultimately, both 1 = 3/2 and 1/2 channels can be expressed by using diagrams A, H, and 
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FIG. 2: Quark propagators corresponding to Eqs. A11-A22 in Appendix A. 
X as follows: 

M(/ = 3/2) - A-X, (10) 
M(J= 1/2) = A-^+ix (11) 

A, H, and X, are schematically shown in Fig. [3] and are given in terms of the quark 
propagators G as 

A(xi,X2,x ll x 2 ) = Ti (G(x 1 ',x 1 )'j 5 G {s) (x 1 ,x 1 ')'j 5 ) x Tr {G(x 2 ', x 2 )j 5 G(x 2 , x 2 )j 5 ) 
= Tr (7 5 G (s) (a;i,a;i')75G'(a;i / ,a;i)) x Tr (7 5 G(x 2 , x 2 ')j 5 G(x 2 ', x 2 )) , 

H(x 1 ',x 2 ',x 1 ,x 2 ) = Tr (G(x 2 , x 2 )j 5 G(x 2 , Xi)7 5 G (s) (zi, Xi) f y s G(x 1 t , x 2 ')^ 5 ) 
= Tr (75G (s) (xi,Xi')75G , (xi / ,X2 / )75G(x 2 ',X2)75G ? (x2,Xi)) , 

X(x 1 ',x 2 ,x 1 ,x 2 ) = Tr (G(xi',X2)7 5 G(x2,X2 / )75G'(x2 / ,Xi)75G' (s) (xi,Xi / )75) 

= Tr { l5 G^(x 1 ,x l ') l5 G(x 1 ',x 2 ) l5 G(x 2 ,x 2 ') l5 G(x 2 ',x 1 )) . (12) 

Here, Tr stands for the trace over color and Dirac indices. 
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H 



X 



FIG. 3: Diagrams A, H, and X. 
C. Calculation of quark propagators using noise vectors 

Now, we calculate the four-point functions in spatial momentum space using the Fourier 
transform of Eq. ( TT2l) at fixed t. In a standard lattice QCD simulation, quark propagators 
G(= D^ 1 ) are calculated by inverting the quark matrix D, 

DX = B, (13) 

using a conjugate gradient type solver. While calculating the hadron four-point functions, 
we obtain the form 

e^Tr [D-\x, «;..)•■■ D~\.., x, t)] . (14) 

X 

For example, in the case of meson-meson scatterings composed of four quark lines, four 
quark propagators, D^ 1 , appear inside Tr[- • •] in Eq. (Jl4"I) . 

If we calculate all necessary components of D~ x with Eq. ( ITBl) . a huge computational 
resource is required. In order to reduce the simulation cost, we introduce noise vectors 

Nr 
3=1 

and rewrite Eq. ( IT41 as 

^ EE e ^ Tr U^D-'frt; ..) • ■■D~\..,y,m j (y)] . (16) 

j x 

Then the four-point functions in the momentum space can be written as 

A, X, H = Y,y ] Z, (17) 



where 

Y or Z = ■■■D- 1 ---£. (18) 

Explicit formulas are given in Appendix [Bj We describe explicitly where the noise vectors 
are inserted. 

In Eq. (fl4"|) . not only x but also the color and Dirac indices are summed up because of 
Tr. Then one may extend Eq. (1151) to include the color and Dirac degrees of freedom, which 
will reduce CPU time further. However, we do not take this in order to keep signal to noise 
ratio at reasonable levels. 



III. NUMERICAL RESULTS 

The lattice simulations are carried out in the quenched approximation at f3 = 2.230 on a 
12 3 x 24 lattice using an improved Iwasaki gauge action. Hopping parameters k uc i — 0.1560, 
0.1580, and 0.1600 and k s = 0.1570 are adopted for these simulations. The lattice spacing 
a obtained using these parameters corresponds to 0.8144 GeV -1 . 

Twenty different configurations separated by 2000 sweeps are used to evaluate the cor- 
relation functions C n , Ck, and Ck-k at each t. We employ a complex Z 2 noise to represent 
the noise vectors in Eq. f|T5|) . The number of noise vectors in this equation is set to be four 
when the color and Dirac indices are fixed. 

One could extend Eq. ( |T5l) to include the color and Dirac indices. In this case, the 
computational time is reduced significantly, but obtained results suffer from large errors. 



The choice here can be considered as a kind of dilution [17]. 

The obtained CV and Ck are shown in Fig. [H Both these correlations are well reproduced 
by one-pole fitted functions for 9 < t < 12, 

Q = Z iC osh( mi (t-N t /2)), (19) 

with i being K or n. The obtained parameters are shown in Table [D 

A. Diagrams A, H, and X 

As shown in Sec. II, only three different diagrams contribute to the correlation functions 
in the present system. Figure [5] shows the results obtained for diagrams A, H, and X for 
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FIG. 4: C n and C K for K q = 0.156, 0.158, and 0.160. 



TABLE I: Parameters for meson propagators 





Zi 


rriia 


CV (Kg = 0.156) 


( 0.35369 ± 0.0002283 ) xlO" 6 


0.722183 ± 0.0004114 


C K {Kq = 0.156) 


( 0.457281 ± 0.0003214 ) xl0~ 6 


0.701035 ± 0.0004315 


C w (Kg = 0.158) 


( 1.08796 ± 0.0008838 ) xl0~ 6 


0.628463 ± 0.0004891 


C K (k q = 0.158) 


( 0.803039 ± 0.0005946 ) xl0~ 6 


0.653207 ± 0.0004654 


C W (Kg = 0.160) 


( 3.85386 ± 0.003642 ) xl0~ 6 


0.526125 ± 0.0005567 


C K (Kg = 0.160) 


( 1.41634 ± 0.001258 ) xlO" 6 


0.606525 ± 0.0005159 



u and d quarks of which hopping parameter, n q = 0.1580. The X and H diagrams become 
negative at large t. At higher values of K q = 0.160, the H diagrams show similar behaviors. 
Since the contributions of H and X diagram distinguish 1=1/2 from / = 3/2, precise 
measurement of these contributions in large t regions is important. 



B. Lattice artifact 



Because the lattice extent is finite, our Kir amplitude contaminated with the artifact 
shown in Fig. [6j These diagram contributions are large in meson correlator calculations at 
finite temperature in lattice QCD, and may mislead us as pointed out in Ref. 18] . Such 
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FIG. 5: Numerical results for A, H, and X contributions in the Kit four-point function at 
K q = 0.158. 

diagram are also seen in the current lattice QCD simulations of the two-meson state. This 
can be easily seen by evaluating the contribution of the fake diagrams in Fig. O 

A ^ e -m„(N t -t) x AKe -m K t + Awe -m„t x ^-m^N^t) 

= 2A n A K e- m * Nt/2 e- mKNt/2 cosh((m K - m n ){t - N t /2)). (20) 

When the mass difference between 7r and K is small, it acts as a constant mode and distorts 
the four-point function at large t. In Fig. [TJ we show the above-mentioned contribution, Eq. 
(I2"0"j) . together with the numerical data corresponding to the four-point functions. 

Another method of avoiding the above mentioned fake diagram is to impose a Dirichlet 
boundary conditions. But here we use a simpler method than it ; we subtract the contri- 
bution of these diagrams numerically from the obtained quantity. A n , Ax, fn^, and nix in 
Eq. (20) can be evaluated from tt and K propagator measurements. In the present case, 
we have the data corresponding to these two-point functions are sufficiently precise to allow 
subtractions of these effects, and hence, Eq. ( 1201 can be subtracted from the Kn four-point 
data. In the following fitting processes, the Kir four-point data is used after the subtraction. 
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FIG. 6: Diagrams that give rise to fake effects. 



C. Fitting analyses 

In general, propagators Cj are composed of many excited states of the same quantum 
number. If t is sufficiently large (1 << t << N t ) and the excited states have significantly 
higher masses than the ground-state mass, the contribution of the lowest energy state is 
predominant in the case of this propagator, and the one-pole model, 



fits the numerical data well. If there are contributions from higher states as well, a two-pole 
model, 



would be more suitable than the one-pole model. Here, mi < rri2, where mi is the lowest 
mass, and m 2 represents contributions from higher states. If the contribution from the 
higher states is very small, the fitting procedure for the two-pole model becomes unstable. 

For K and n two-point functions, the one-pole ansatz, Eq. ()2ip . works well, and hence,the 
masses of K and tt are obtained with high accuracy (Fig. Hj). The results of our simulation 
have already been shown in Table [H and the values are consistent with those provided by 



For calculating the propagators of the Kn system (meson four-point function), we adopt 



C;(0, t) = Z 1 cosh(-mi(t - N t /2)) 



(21) 



Ci(0, t) = Z 1 cosh(-mi(t - N t /2)) + Z 2 cosh(-m 2 (t - N t /2)), 



(22) 



CP-PACS. 
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with lattice artifact 
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1= 1/2 

without lattice artifact 
1 = 3/2 

with lattice artifact 
1 = 3/2 

without lattice artifact 




1=1/2 

with latice artifact 
1= 1/2 

without lattice artifact 
I = 3/2 

with latice artifact 
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FIG. 7: Kir four-point functions before and after subtracting the lattice artifact, Eq. (|20[) . The 
figures from top to bottom correspond to K q values of 0.1560, 0.1580, and 0.1600, respectively. 
In each figure, for both / = 1/2 and 3/2 channels, naive four-body correlations, Kn, and those 
obtained after correction for the lattice artefact, Eq. (|20p . are shown. 

a two-pole model by taking into account higher excited states; however, the calculation 
method in this case is not simple. A naive application of Eq. ( |22i) is not stable when the 
fitting region we use is changed. In order to obtain reliable results, we take the steps as 
follows: 

1. We apply Eq. ( 1221) to Ck w by changing the fitting region. 

2. We choose a stable region, where the obtained mi shows a plateau, and calculate the 
average of mi in this region. 

3. We also apply one-pole fitting and verify that mi obtained by two-pole fitting is lower 
than that obtained by one-pole fitting. 
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4. If the fitting procedure for the two-pole model is unstable, we adopt the results ob- 
tained with the one-pole model. 

In the present simulation, N t is 24 and the source field is set at t — 0. Because of its 
bosonic property, the propagator is symmetric at t = N t /2 = 12. We apply the two-pole fit 
to propagator Ck-k in the region t s < t < N t /2. The larger t s corresponds to a propagator 
in the larger t range, which is more reliable for picking up the ground state. However, in 
this case the number of available lattice points decreases. Figure M shows m\ as a function 
of t s . Then, we obtain the statistical average of the plots in Fig. [S]in order to obtain our 
final results, which are shown as horizontal lines in the figure. 



>ne-pole fit for region 9 to 12 




- • k=0.158 1=1/2 



FIG. 8: Effective mass extracted by one-pole and two-pole fittings as a function of ts- 1 = 1/2 
and k = 0.158. The central region (4 < t s < 8) is used to determine the ground-state energy. The 
value obtained and its error are shown by horizontal solid and dashed lines, respectively. The data 
point on the vertical axis shows the result obtained by the one-pole fit. 

For I = 3/2 with k = 0.1560 and 0.1580, two-pole fitting yields results with large statis- 
tical errors. Therefore, in these cases, we adopt one-pole fitting in the region 9 < t < 12. In 
the other cases, two-pole fitting gives better results than one-pole fitting, i.e., mi is smaller 
than the mass obtained with the one-pole model, and the statistical error is sufficiently 
small. The obtained values of Ek-k and AE = Ek-k — ("V + mic) at different values of the 
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TABLE II: Values obtained by fitting procedures for four-point function Ckk and two-point func- 
tions Ck and C n . The S-wave scattering length ao is also calculated using these values. 





^"ud 


Ek-tt 


5Ex-k 


AE 


a 




0.1560 


0.6542 


0.01804 


-0.7690 


0.4532 


I = 1/2 


0.1580 


0.9264 


0.01882 


-0.3552 


0.1598 




0.1600 


1.066 


0.008285 


-0.0672 


0.0222 




0.1560 


1.379 a 


0.02503 


-0.0446 


0.01850 


I = 3/2 


0.1580 


1.258 a 


0.03848 


-0.0240 


0.00896 




0.1600 


1.087 


0.002129 


-0.0453 


0.01505 



a One-pole fit 



hopping parameter k are summarized in Table [TT1 
D. Chiral extrapolations and scattering length 

In the present study, the pion masses m n are considerably larger than those determined 
experimentally, and hence, we need to adopt an extrapolation procedure. Although the 
behavior of E] <n near the chiral limit is not very clear, we plot E 2 Kv for I = 1/2 and 3/2 as 
a function of 1/n in Fig. [TOl together with m 2 , and m 2 K ; this is because Exn is expected to 
be dominated by m n and tuk- 

The least mean square procedure for a linear function, f(x) = ax + b, provides us the 
errors for both a and b. The least square procedure assumes that the fitted line passes 
through the center of the weight of plotted points. Therefore, errors of a and b are not 
independent but they are correlated; in the case that a fluctuates positively b must fluctuates 
negatively and vice versa. Hence, we may denote the line y = ax + b as the middle line, y = 
(a + Aa)x + (b — Ab) as the lower line and y = (a — Aa)x + (b+Ab) as the upper line, with Aa 
and Ab being the error of a and b, respectively. With these three lines, we evaluate the errors 
in the extrapolation process. Putting physical pion mass as = 139.57 MeV, the fitted 
line of m 2 provides us hopping parameter at physical point as 1/^^=6.077724 ±0.003071. 
Here, the error, A(l/n phys ) = 0.003071, is one half of the difference between the value of 
the upper line and the one of the lower line at the physical pion mass. Also in the case of 
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* C Kii (k=0.156, 1=1/2)* 10 
fitted curve 



C Kji (k=0.158, 1=1/2) 
■ fitted curve 



▼ C Kji (k=0.160, 1=1/2) /10 
fitted curve 



C K7t (K=0.156, 1=3/2) * 10 
■ fited curve 

C Kii (k=0. 158, 1=3/2) 
- fitted curve 



▼ C Kjt (K=0.160, 1=3/2)/ 10 
fitted curve 



FIG. 9: Four-body correlator (propagators of the Kir system ) and fitting curves for different k 
values. Upper and lower figures correspond to the 1 = 1/2 and / = 3/2 channels, respectively. 

the energy of the two-particle state, we denote three curves as the middle, the upper and 
the lower which are the square root of the corresponding lines of the linear fit to the square 
of the energy of the two-particle state, respectively. The central value at physical points is 
determined by the middle curve and 1/K p h ys - The error originates from the linear fitting, 
Ai, is estimated from the difference between the upper curve and the lower curve at 1/K p h ys . 
There exists another kind of the error, A 2 , which originates from 1/ ' n p h ys which is evaluated 
from the difference of the middle curve at 1 / n p h ys + A ( 1 / n p h ys ) and at 1 / n p h ys — A ( 1 / K p h ys ) • 
As a final error value, we may take root mean square of two independent errors, Ai and A 2 , 
as AE(rc,K) = yAf-P~Af. We also estimate errors of E(ir,K) in the other channel and 
ttlk in the same manner. 

Substituting the obtained values into Liischer's formula, Eq. ([T|), we can obtain scattering 
length ao- Unfortunately, the results obtained in the present study are limited to a fixed 
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6 ~/6A 6.2 6.3 6.4 6.5 
1/Kr 1/k 



• (m n a) 

o (m K a) 2 

■ (E K7l a) 2 : 1=3/2 

(E Kjr a) 2 : 1=1/2 



FIG. 10: Chiral extrapolations of m 2 , rn 2 K , and E\ . The horizontal axis is 1/k,j, where K q is the 
hopping parameter for u and <i quarks. The vertical line represents 1/k c = 6.069351464. 




6 /f l 

1/Kc Physical Point 



6.2 6.3 6.4 6.5 

1/K 



▼ m K a 

■ AE KjI a : 1=1/2 
♦ AE KjI a : 1=3/2 



FIG. 11: AE = Exn — ( m n + thk) as a function of 1/k 9 . 

volume; nevertheless, by using the physical size of the lattice unit, we can evaluate do- 
Let us rewrite Eq. (OQ) as 

(23) 
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TABLE III: qq at the physical point where m 7r =139.75 MeV. 







a (MeV) 


a m n 


I = 
I = 


1/2 
3/2 


-0.8794 

fl 1 1 7R +0-0712 
-U.ll/8 _ .0901 


-0.6248 +8:ggi 

n +0.0506 

-U.U86( _ .0640 



where 

^(^)=^[l + cif + c 2 (^) 2 ]+0(L-). (24) 

The left-hand side of Eq. (|23l) represents the shift in energy caused by the interaction between 
7r and K; this energy shift can be evaluated from the correlations. The right-hand side, 
f2 fsp.) , represents the effect of a in the unit of box length. Figure [T21 shows fl as a function 
of ao/L. Q changes slowly at around ao/L ~ 0, indicating that a small ao is sensitive to 
AE at around Ai? ~ 0. Further, ao easily changes from a positive small value to a negative 
small value and vice versa. 

The chiral extrapolation of ao is shown in Fig. [T3j Our final results are summarized in 
Table IIHI At all the calculation points, a is positive, as shown in Table [III however, the 
chiral extrapolation for E\ v scaling leads to a change in the sign of AE and a subsequent 
change in the sign of ao- At the physical point where m n = 139.75 MeV, ao(I = 3/2) 
= —0.0837 and a^{I = 1/2) m n = —0.6248. In both cases, AE > 0, the phase shift 
5 < 0, and forces are repulsive. 




FIG. 12: £1 evaluated through AE obtained by fitted curve in Fig. [TT] (left-hand side) and Q as a 
function of ao/L (right-hand side). 
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FIG. 13: S-wave scattering length ao as a function of 1/K q . Simulation data points at three values 
of Kg and extrapolated values at = 140 MeV. 

IV. CONCLUDING REMARKS 



The scattering length, qq of the I = 3/2 channel in scatterin g ha s been studied 



by theoretical and experimental approaches so far 13, |2Cj, |2l|, |22|, |23|, |24j, |25| . Previous 
experiments have reported that a^^m^ has a small negative value, i.e., —0.13 ~ —0.05 19, 



20 



21] . Small negative value was also claimed by a theoretical model based on the chiral 



perturbation theory as a = —0.129 



-0.05 
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The first lattice calculation of Kir 



scattering in the 1 = 3/2 channel was performed by Miao et al. [lOj . and the value of ao was 
found to be -0.048. 

On the other hand, no direct simulations have been carried out on Kir scattering in the 
1 = 1/2 channel for estimations of ao- The NPLQCD group carried out lattice simulations 
for Kit scattering in the 1 = 3/2 channel and evaluated the ao values for both 1 = 3/2 
and 1/2; the results of / = 1/2 were obtained on the basis of chiral perturbations and not 
directly by simulations 11]. They obtained a small negative value of —0.0574 for the 1 = 3/2 
channel and found that = +0.1725^ l ' )l j ) 57 for the 1=1/2 channel. The pion mass mo- 

used by them was lower than that used in our present study. We could carry out a direct 
comparison of our results with those of NPLQCD group if we carry out our simulations at 
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a low m^. Flynn and Nieves|26|] used the scalar form factors in semi-leptonic pseudo-scalar- 
to-pseudo-scalar decays to extract information on Kit scattering in the 1 = 1/2 channel and 
obtained m n a 1/2 = +0.179(17) (14). 

In this paper, we presented a lattice QCD simulation of the Kit scattering length and 
formulas that make this calculation simple and feasible. We have observed the followings: 

• I = 1/2 and 3/2 amplitudes can be expressed using quark diagrams, A, H and X. 

• use of a reasonable noise vector works well. 

• a simple method for eliminating fake diagrams caused by the finite size of lattice is 
useful. 

• the accurate measurements in large t regions is important. 

• the behavior of AE in the chiral extrapolation needs careful analysis. 

Thus, it is now possible to study Kir scattering reactions on the basis of lattice QCD 
simulations. However, our present study does not provide adequate information on scalar 
mesons with the strangeness, i.e., k, because our study is restricted to the zero- momentum 
case. Therefore, this paper does not include results on the phase shift 5{p), which is indicative 
of a k pole in the Kir channel. The study of Kit interactions is the first step in the study of 
hadron interactions including s-quarks, and we are now entering into the new era to study 
hyperon interaction from QCD. 

We show exchanged mesons for meson-meson and nucleon-hyperon interactions in Table 
IIVI and [V] 28j NN interactions include 7r exchange and require lattices of very large 
size for estimation of the its large scattering length. On the contrary, NA interactions 
can be studied using a lattice of reasonable size. Hence, NA interactions are suitable for 
lattice studies because they have interaction ranges that can be fit to the Liischer's formula; 
these interactions will be extensively studied in future experiments in J-PARC and other 
laboratories. 
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TABLE IV: Exchanged pseudo-scalar mesons for baryon-baryon scatterings. 



scattering 


system 


exchanged mesons 


N - N 




TT, 7], rf 


A-N 




K, 77, rf 


E - N 




TT, K 


E-N 




TT, 7], r( 



TABLE V: Main contributions in t/u- and s-channels for pseudo-scalar meson- meson scattering. 



scattering system 


t/n-channel 


s-channel 


TT — TT 
K-TT 


p, a 
p, a 


p, a 

K*, K 
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APPENDIX A: I = 1/2 Kir SCATTERING AMPLITUDES IN TERMS OF 
QUARK PROPAGATORS 



(OKo(V)o, + (x 2 ooyzi)o; + (x 2 )) 

= Tr (G^(x 2 ,x' 2 ) l5 G^ (Sarahs) Tr (G^^x^G^x^x',)^) (Al) 

- Tr (GW(x 2 ,x 1 ) 75 G' (s H^^i)75G' (d) (4,4)75G (u) (x 2 ,x 2 )75) (A2) 

(O x o(x 1 ')O ff+ (a; 2 , )0^ + (a;i)Oio(x 2 )) 

= --^Tr (GW(i 2i i 3 ) 75 ) Tr (G^^, ^G^i, 4h 5 G (<i) (4 4) 7s) (A3) 

+^Tr(GM(x , 2 ,x 2 )75G (tt Hx 2 ,x 1 )75G (a H^i^075G w (a/ i y 2 ) 7 5) (A4) 

+ _L T r (GW(x 2 ,x 2 ) 75 ) Tr (G«(a;i, 4)75^^(4, 4)75^^(4, xi) 75 ) (A5) 

_^Tr(GW(x;,x 2 ) 75 GW(a: 2 ,x 2 )75G("H4^i)75G (s) (x 1 ,x / 1 )75) (A6) 

(O jr+ (a; 1 , )0 w o(a; 2 , )0^,(xi)0i + (x 2 )) 

= ~Tr (G^(4,4) 75 ) Tr (G (u) (:4 x 2 )7 5 G (d) (x 2 , Xi)7 5 G (s) ( a: i> 4)75) (A7) 

+^Tr(GW(x , 2 ,x 2 )75GW(x 2 ,xi)7 5 G( 8 H^i,4)75G( u H4,4)7^ (A8) 

+ _L T r (GW(x' 2 ,x' 2 ) 75 ) Tr (G (d) (x 2 ,x 1 )75G w (^i,4)75G (u) (x , 1 ,x 2 )75) (A9) 

-^Tr(GW(4,x 1 ) 75 GW(x 1 ,x 1 075G( tt H4^2)75GW(x 2 ,x , 2 )75) (A10) 



23 



(O K+ ( Xl ')OAx2')O^ + (xi)Ol (x 2 )) 

04 4) 75 ) Tr (GW(x 2) x 2 ) 75 ) Tr (G^(xf 1 ,x 1 ) l5 G^(x 1 ,af 1 ) l5 ) (All) 

(x' 1 ,x 1 h 5 G^(x 1 ,x[)^ 5 ) Tr (GM(x' 2 ,x 2 )kGW(x 2 ,x' 2 )k) (A12) 

(x' 1 ,x^)75G^(4^2)75G {M) (x 2 , a : 1 )75G (s) (^i^ / i)75) (A13) 

(x / 1 ,4)7 5 G^(4^i)75G (s) (x 1 , a ;' 1 )75) Tr (G^(x 2 , x 2 ) lb ) (A14) 

K,x 2 )7 5 G' (M H^,4)75G' {M) (4^i)75G' (s) ^i^ / i)75) (A15) 

(x / 1 ,x 2 )7 5 G^(x 2 ,x 1 )7 5 G' (s) (a; 1 , a :' 1 )75) Tr (CW(4, z 2 ) 75 ) (A16) 

■' 1 ,x 1 h 5 G^(x 1 ,x' 1 )^ 5 ) Tr {G {u \x' 2 ,x' 2 ) l5 ) Tr (G«(i 2) z 2 ) 75 ) (A17) 

(x[,x' 2 h b G^(x' 2 ,x 1 h 5 G^(x 1 ,x' 1 ) l5 ) Tr (G (<f) (i 2 ,x 2 ) 75 ) (A18) 

(x' 1 ,x 1 h 5 G^(x 1 ,x > 1 )^ 5 ) Tr (G^(x 2 , x 2 ) 75 ) Tr (G^(x' 2 , x' 2 ) lb ) (A19) 

(x' 1 ,x 2 ) 7 5G {M) (a; 2 ,x 1 )75G' {s) ^i^ , i)75) Tr (G (d \x' 2 , x' 2 ) l5 ) (A20) 

(4,4)75) Tr {gW(x 2 ,x 2 )k) Tr (G^i, ^hsG^Vi, xi) 75 ) (A21) 

(x' 2 ,a: 2 )75G (d) (^2,4)75) Tr (G^, x[) l5 G^ (x[, Xl ) l5 ) (A22) 



-^Tr 

+^Tr (G (M)i 
-^Tr 
+^Tr (G^ 
-^Tr (G {u 
(G {u: 

2 



Hi 



-^Tr (G (u 
+^Tr (G (M) i 
-^Tr 

+^Tr (G& 



APPENDIX B: EXPLICIT FORMS OF THE DIAGRAM A, H AND X 

In this section, we employ two independent random noises, £ and 77. They are a function 
of the spatial coordinate, i.e., they live on a each time slice. 



N R . 

— VjixYvM = <W 



(Bl) 
(B2) 
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1. Diagram A 



A (p 

rii VniVmi Pm ) (B3) 
. . „+ip„-xi „-ip n -x 2 ' „+ip m -xi p-ipm-32 



x A((x 1 ',t),(x 2 ',t),(x 1 ,t s ),(x 2 ,t s )) (B5) 

-ipn-xi' e +ip 

T 3 W 



p + iPn-Xl' p+iPm-Xl 



x\' x\ 

x Tr ( 75 G (s) (fi, t s -x^t) lh G(x^t-x x M)) (B7) 

D -ip n -x 2 ' e -ip m -x 2 

~T 3 IX 



x Tr(7 5 G , (f 2 ,t 5 ;^',t)75G'(f 2 ',t;f 2 ,t5)) (B9) 
ii j/i xi' xi 

x Tr (4(y 1 ) 75 GW(y 1 ,t 5 ;5 ? i , ,t)75G'(fi , ,t;fi,t s )0 1 (5 ? i)) (Bll) 

x J2J2J: e —— < B12 > 

32 V2 X 2 ' X2 

x Tr (vli&hsGifa, t s ;S 2 ',th 5 G(x 2 ',t;x 2 ,t s ) Vh (x 2 )) (B13) 
We can write the trace terms as 
Tr ( 75 G' (s) (fi, t s ; xi, t) 75 G(fi', t; x u t s )) = ^(a, a\-f 5 G (s) {x u t s ; xi, t)^ 5 G{x /, i; f i, t s )|a, a) 

a, a 

(B14) 

Tr (7 5 G(f 2 , *s; 3/, t)^ 5 G(x 2 , t; x 2 , t s )) = y~]{a, a\~f 5 G(x 2 , t s ; x 2 , t)~f 5 G(x 2 , t; x 2 , t s )\a, a) 

(B15) 

where a and a stand for the color and Dirac indices, respectively. 
Then 

1 y _ e +ip n -xi' ^ 

A(p n ,-p n ,p m , -p m ;t,t s ) = J^z2 z2z2 z ^^ t '^ t sy — 7-3 — Yi(xi,t',t s ) 



ji a,a xi' 

1 p—lpn-X2' 

EEw*;^- fjr-YiWMs) (Bi6) 



j2 a,a x 2 ' 
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where 

Y 1 = G— or (B17) 

Y 1 (x 1 ',t;t s ) = ^Gfa'^Zuts) ^ ( B18 ) 

Z x = l5 G^^ h \a,a) or (B19) 

Z^x^t-ts) = ^7 5 G' (s)t (fi , ) t;yi,t5)75^ 1 (yi)|a,a) (B20) 

Y 2 = G^—rj j2 \a : a) or (B21) 

n(x 2 ',t;*s) = ^G(x 2 ',t;x 2 ,ts) ^ 77 j2 (x 2 )|a,a) (B22) 

Z 2 = 7 5 G' t 7 5 r7 j2 |a,a) or (B23) 

Z 2 (x 2 ',t;t s ) = J2^^2,t;y 2 ,tsh 5 V 3 M)\a,a) (B24) 

m 

where f or Hermite conjugate includes the color, Dirac and site indices. 
Indices of these vectors are 

Y^,t;t s ) =X;X;G^(f,*;f,* s )e(5)|a,a) C)7 (B25) 

x c,7 

Zj(f ,t;*s) = (0 5f ,t s ;f ,i))*^)k,a) C;7 (B26) 

^ • f = E E E *; * *; t s ) (B27) 

b P x' 
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2. Diagram H 



H (Pn,-Pn,Pm,-Pr, 



EEEE 



g-ip-n-Sl' e +iPm-Xl g"«Pm-^2 



L 3 L 3 L 3 L 3 

Xi' X2' Xl X-2 

x H((xi, t), (x 2 , t), (x^ts), (x 2 , t s )) 



L 3 L 3 L 3 L 3 



x Tr (75G (s) (x 1? i s ; xi, t)-f 5 G(xi, t; x 2 , t)-f 5 G(x 2 , t; x 2 , t s )lsG(x 2 , t s ; x u t s )) 

! e +iPnSl' e -iPnS2' e +ip m Sl e -ip m S 2 

l^l^l^l^l^ us l 3 L 3 L 3 

j yi xi' x 2 ' xi x 2 



x Tr [^(mh^^yuts^x^t^Gix^t^x^t) 

x 75G(f2 , ,t;f 2 ,t5)75G'(f 2 ,t 5 ;fi,t 5 )^(fi)) (B28) 



We can write a trace term as 



Tr (7 5 G (s) (fi, t s ; x x ' , t)i 5 G(xi, t; x 2 , t)~f 5 G(x 2 , t; x 2 , t s )lbG{x 2 , t s ; x u t s )) 
= y^(a, a|7 5 G (s) (fi, t s ; xi, t)-f 5 G(xi, t; x 2 , t)~f 5 G(x 2 , f; x 2 , t s )^ 5 G(x 2 , t s ; x u t s )\a, aQB29) 

Then 

1 . _ g-ipn-Xl' _ 

H(Pn,-Pn,P m ,-Pm]t,t S ) = z2 ^ , t] t S ) j£ Y(x 2 ',t;t S ) (B30) 

R j Si' 

where 

e ~iPm-X2 p+iPm-Xl 

Y = L3— ^3-^1°' »> or ( B31 ) 

^ e +iPm-Xl 

Y(x 2 ',t;t s ) = 22^2G(x 2 ',t;x 2 ,t s ) — ^ — lsG(x 2 ,t s ;xi,ts) — ^ — 

X2 Xl 

x ^-(fi)|a,Q!) 

= J2 G (^',t;x 2 ,tsh 5 ^ yi(£ 2 , *g;te) (B32) 

£2 

Here 

L 3 



Yi(x 2 ,t s ;t s ) = 2jG(x 2 ,tsf;fi,*s) — 7^ — ^(fi)|a,a) (B33) 
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Z = ^7 5 G f — ^75G (s)T 7 5 0|a,«) or (B34) 

Z(x 2 ',t;t s ) = 7sg t (^ / , t; f 1 / , t) ^ ^ 75 G (a) t; y u t 5 ) 75 (B35) 
x &(j/i)|a,a) (B36) 

p-iPn-Sl' _ 

= Y^lsG^x^t-x^t)— — Z^x^t-tg) (B37) 

Xl' 

Here 

fi(fi',*;*s) = ^T5G (a)t (xi',t; j7i,ts)75^(j/i)|a,a) (B38) 
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3. Diagram X 



X (Pn,-Pn,Pm,-Pm) 

e +iPnSl' e ~ipn-X2' e + ip m -Xl e ~ip m -X2 

1^1^ 1^1^ L 3 J} J} L 3 

xi ' a?2' xi X2 

X -X"((^l', *), fa', t), (f i, f S ), (£ 2 , 

= EEEE 



e +ip„'X*i e ~ip n -X2 e +ip m -xi e ~ip m -X2 



L 3 L 3 L 3 L 3 



X 



Tr (75G {s) (fi, t s ; £/, t)7 5 G(fi / , t; £ 2 , ts)l*>G{x 2 , t s \ x 2 , t)^ 5 G(x 2 ', t; x u t s )) 

1 + ipn-X\' — ip n -X2' +ipm-Xl c —ipm-X2 

l^l^l^l^l^ l 3 T 3 T 3 L 3 

j J/l Xl' X2 1 XI X2 

x Tr (eJ(yi)7 5 G' (s Hfi,^;^ , ,t)75G'(fi , ,t;f 2 ,t5)75G(f 2 ,t5;x2',t) 

x 7 5 G(f 2 , ,t;fi,t 5 )0(^i)) (B39) 

We can write a trace term as 

Tr (7 5 G (s) (fi, t s ; t)i 5 G(xi, i; f 2 , t s )^G{x 2 , t s ; x 2 , t)~j b G{x 2 , t; fi, t s )) 
= y^(a, a|7 5 G (s) (f i, i s ; f /, t)n/ 5 G(x i, t; f 2 , t s )^ 5 G(x 2 , t s ; x 2 , t)~f 5 G(x 2 , t; x u t s )\a, a£B40) 



Thus 



*(Pn, -Pn, An, "Pm) = ^" ^ S Z ^ 2 '' t; *s) f — p~ ^l^', *s) (B41) 

^ J a,a X2 1 
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Here 

Y 1 (x 2 ',t;t s ) = 2_,G{x 2 ',t;x u t s ) ^ Zj{xi)\a,a) (B42) 
Z(x 2 ',t;t s ) = 2_,2^2^-i b G ] (x 2 \t]x 2) ts) — — z — ^G\x 2 ,t s ;xi,t) 

x 2 Si' yi 



x — ^— 75G W W, f; ft, ts) 7 50(yi)|a, «> 



^ 75^(^2', *; x 2 , -— f '(f 2 , * s; i) (B43) 



L 3 

Z'(x 2 ,t S ;t) EE J] J] 75 Gt(f 2 ,t s ;f 1 ', f )_^_ 75 GW T (f 1 ',t;j7 1 , is ) 75 
xi ' j/i 

x &(yi)ka) 

g-ipn'Xl' ^ 

= 2^ 75^(^2, *s; fi', f) — ^ — 1', t; t s ) 

XI ' 

= ^75G wt (fi',*;j/i,ts)75^(j/i)|a,a> (B44) 
2/1 
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